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We simulate a phase-separating bilayer in which the leaflets experience a direct coupling favour¬ 
ing local compositional symmetry (“registered” bilayer phases), and an indirect coupling due to 
hydrophobic mismatch that favours strong local asymmetry (“antiregistered” bilayer phases). For 
wide ranges of overall leaflet compositions, multiple competing states are possible. For estimated 
physical parameters, a quenched bilayer may first evolve toward a metastable state more asymmetric 
than if the leaflets were uncorrelated; subsequently, it must nucleate to reach its equilibrium, more 
symmetric, state. These phase-transition kinetics exhibit characteristic signatures through which 
fundamental and opposing inter-leaflet interactions may be probed. We emphasise how bilayer phase 
diagrams with a separate axis for each leaflet can account for overall and local symmetry/asymmetry, 
and capture a range of observations in the experiment and simulation literature. 


INTRODUCTION 


Phase separation within mixed bilayers is intimately 
linked to the symmetry or otherwise of their separate, yet 
coupled, leaflets |ll-[llj. The overall compositions of the 
leaflets may by asymmetric; for instance, one leaflet con¬ 
taining an equal mixture of saturated and unsaturated 
lipids and the other predominantly unsaturated lipids. 
The phase behaviour of such bilayers differs dramatically 
from those whose overall leaflet compositions are sym¬ 
metric [H . Alternatively, even if the overall compositions 
are symmetric, as is common in model systems, phase 
separation can lead to either registered or antiregistered 
phases, and thus local symmetry or asymmetry [ 12 . [l^. 

There is evidence that bilayers experience competing 
inter-leaflet coupling effects. Many observations have 
shown separation into registered (R) bilayer phases, those 
which are locally symmetric, comprising leaflets with the 
same phase and composition |2l44l.ll3l|. This implies direct 
inter-leaflet composition coupling via a mismatch energy 
per area 7 S i S El El- Such a direct coupling is 
also used to explain apparent induction or suppression 
of domains in one leaflet by the other, in cases where 
overall leaflet compositions are asymmetric [ 2 , 
timates from theory [Bl, @, simulation El El 
more recently, experiment vary over approximately 
7 = 0.01 — l/cB 2 "nm“^. Conversely, lipid hydrophobic 
length mismatch indirectly couples the leaflets via their 
combined hydrophobic thickness. This effect promotes 
strong local asymmetry, i.e., antiregistered (AR) bilayer 
phases, to give more uniform bilayer thickness and avoid 
the elastic cost El that would otherwise result from hy¬ 
drophobic mismatch [ 1201 ^. | 2 Q| - [^ . 

A microscopic model [y incorporating these effects pre¬ 
dicts a kinetic competition of R and AR bilayer phases. 
Linear stability analysis in physical parameter ranges 
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FIG. 1. (Color online) Schematic free-energy landscape 
f{(jf,(j)^), which determines available phase coexistences [i|. 
Cartoons show the dominant transbilayer arrangement of 
model S and U species in R and AR bilayer phases. Curved 
arrows represent competing linear instabilities for a bilayer of 
(0.5, 0.5) overall composition. 


predicts that instabilities to forming registered or an¬ 
tiregistered phases can be of comparable strength, so that 
which occurs first after a quench is determined by molec¬ 
ular properties such as tail length or unsaturation Q. 
If a metastable state forms first, nucleation is required 
to subsequently initiate the equilibrium (typically regis¬ 
tered) phases 

In this paper we examine the resulting kinetics by di¬ 
rect simulation of the microscopic lattice model [l| . Mul¬ 
tiple phase coexistences involving registered and antireg¬ 
istered phases compete over wide regions of the phase 
diagram. Even in cases where the given overall leaflet 
compositions prevent the attainment of either “perfect” 
local symmetry or asymmetry throughout the bilayer, 
the bilayer may still select between states with greater 
or lesser amounts of registered and antiregistered phases 
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than the uncorrelated case of fully-independent leaflets. 

We find kinetic signatures of an Ostwald stage rule 
by which a bilayer progresses through metastable states 
on its way to equilibrium. For example, the nonmono¬ 
tonic evolution of a parameter describing the degree of 
molecular transbilayer symmetry could be measured in 
molecular simulation or in experiment. We emphasise 
the utility of bilayer phase digrams with a separate com¬ 
position axis for each leaflet [1, [H, [H . We show how these 
relate to available experimental results, and to conven¬ 
tional phase diagrams that do not account for transbi¬ 
layer symmetry / asymmetry. 


II. LATTICE MODEL 

The model introduced in [l| describes a local bilayer 
patch as N sites on a square lattice of spacing a ^ 0.8 nm, 
where each site contains top and bottom leaflet lipids. 
Each lipid has a hydrophobic length from which we 
define the bilayer thickness 

di=£\+£^, ( 1 ) 

and leaflet thickness difference 

\=£l-£^. ( 2 ) 

Model species S and U represent either saturated and 
unsaturated lipids, or the more ordered (liquid-ordered 
Lo or gel) versus liquid-disordered {Ld) states of a ternary 
mixture f Section III Pp . 

Defining = 1 if the top (bottom) of site i contains 

an S lipid, = 0 if I/, the Hamiltonian is 

^ = E + E 

<i,j> ^ ' <i,j> 

+ ^ ^ iK {{£\ - + {£\ - , 

i i 

(3) 

where species-dependent ideal (i.e., preferred) hydropho¬ 
bic lengths are ^so for an S lipid at the top 

(bottom) of site i, or £uo for U [25|. Each site is pair¬ 
wise registered {SS or UU) or antiregistered {SU oiUS). 
An R bilayer phase is one dominated by SS or UU sites 
such that, microscopically, most lipids face one of the 
same type in the apposed leaflet, while in an AR phase 
most lipids face one of the opposite type. 

The Ising or Elory-like Vuv captures length- 
independent interactions occurring among neighbours 
within each leaflet, arising from, e.g., headgroup interac¬ 
tions. Its strength is quantified by E = Eio —|(Eoo+hii)- 
In mean-field, phase separation of the Ising model (i.e., if 
no other interactions were present) requires V > 0.5 /cbF, 
while V > 0.88 /cbF is required in the presence of fluctu¬ 
ations, as in a simulation [26|. 


The “direct” inter-leaflet coupling B is suggested by 
observations of registered domains in experiment and 
simulation lii, El m, as well as the apparent “in¬ 
duction” of domains in one leaflet by those in the other 
[T^ (see Appendix These observations imply a di¬ 
rect composition coupling between the leaflets, resulting 
in a mismatch e nergy per area for local compositional 
asymmetry H, H, El El El El- The direct coimling is 
expected to depend on tail structural features [51,0; here, 
tail structure is implicitly mapped to tail length , and 
the B term in our model promotes pairwise R by penalis¬ 
ing tail length mismatch across the midplane. The exact 
mechanisms (e.g., transmidplane tail interactions, leaflet 
curvature) are not crucial, though, and for comparison 
with the literature we can estimate an effective conven¬ 
tional inter-leaflet mismatch energy 7 (shown later on an 

axis of Eig. [4j) [El, 0, [13, [H, [13, [3 • 

The hydrophobic “indirect” coupling J promotes pair¬ 
wise AR, by penalising mismatches in the bilayer thick¬ 
ness profile [i3,[i3,0- We define J = 4J, which appears 
in the mean-field approximation of Eq. [3] from which the 
coarse-grained free-energy density (Eig. [T|) is derived [l|. 
K can be related to the area compression modulus sa 
by mapping length changes to area changes assuming 
constant volume of the hydrophobic tails (described in 

IH). It quantifies the penalty for deviation from ^ 
so that smaller values soften both inter-leaflet couplings. 
While the interplay of J and s gives a reasonable en¬ 
ergy scale for the line tension associated with hydropho¬ 
bic mismatch it does not capture the detail of the 
bilayer’s internal elastic deformations 19|. We will ex¬ 
plore implications of the findings of [ly in an upcoming 
Comment on that work [ 2 ^. 

Eor simplicity the parameters H, J, and n are assumed 
to be the same for both lipid species. This will not af¬ 
fect the key qualitative features that are of interest here, 
but relaxing this assumption could lead to interesting be¬ 
haviour in specific systems. 

Ao = £so — £uo quantifies the intrinsic mismatch in 
both preferred tail length and structure between species; 
although cast as a tail length difference, it works with 
both J and B to control the strength of the indirect and 
direct inter-leaflet couplings. We model varying the tail 
hydrophobic length mismatch alone by varying J, while 
varying B alone corresponds to changing the degree of 
tail structure mismatch (for instance by increasing the 
difference in unsaturation). A reference total thickness 
do = iso + iuo may be defined, but in the absence of an 
external field acting on bilayer thickness the absolute val¬ 
ues of iso and iuo are irrelevant and only their difference 
plays a role. 

Recent work employed a model capturing similar com¬ 
peting inter-leaflet coupling effects to those studied here 
|23| . The authors successfully analysed local correlations 
and clustering via a quasi-chemical approximation, in the 
context of a laterally homogeneous bilayer rather than 
one undergoing phase separation. 
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A. Mean-field free energy landscape 

Eq. [3] leads to the mean-field local free-energy density 
f{4>\ </>^ d, A) Q as a function of average leaflet com¬ 
positions in the local patch 

A/'t(b) 

^t(b) ^ ^ , (4) 

where there are top (bottom)-leaflet S lipids in the 

local patch, and of thickness variables 

d= , 

( 5 ) 

where the sums are over the local patch. The local free- 
energy density as a function of local composition is 

j[ann.](^t^^b) = ^b^ J[ann.]^ ^[ann.]^ ^ 

The label [ann.] (annealed) indicates equilibration of the 
local thickness variables at given local compositions [l| 
and will be omitted hereafter. Each leaflet is described 
by a separate composition variable, and bilayer phase 
behaviour can then be calculated in the space of ( 0 ^, 0 ^) 
fSection III Bp. similarIv to Hi- An illustrative free- 
energy landscape is shown in Eig. [TJ 

The model species have the same molecular area [l|. 
Breaking this assumption, or otherwise breaking the sym¬ 
metry of the free-energy landscape (e.g., allowing param¬ 
eters such as hi to vary between species) should not qual¬ 
itatively affect the results. (However, in experiment or 
molecular simulation one would typically interpret 
as area fractions, which are equivalent to number frac¬ 
tions in our model. Eor instance, an “equimolar mixture” 
corresponds to an “equal area fractions mixture” in a real 
system.) Such a change would remove the symmetry of 
the free energy (Eig. [T]) and phase diagram (Eig. [ 2 ) 3 , c) 
under reffection through the = 1 — line. A separate 
consideration is that biological membranes with, e.g., dif¬ 
ferent solvent environments for each leaffet could break 
the top-bottom symmetry, thus removing the symmetry 
under reffection through the (j)^ = (j)^ line. 

B. Bilayer phase diagrams and coexistence 

Eor the simulations presented below, the bilayer size is 
£2 _ 200 ^ = M. Script letters refer to the entire bilayer, 
while N used above described a local, homogeneous bi¬ 
layer patch. Similarly, we explicitly define the overall 
leaffet compositions 

iydCb) 

^t(b) _ _ S _ 

AT ’ 

as distinct from the local compositions defined in Eq. 01 
(although the distinction is conventionally clear from 
context M)- 


Coexisting bilayer phases correspond to points on the 
free-energy surface ^0^,0^) (Eig. [T]) that share a com¬ 
mon tangent plane [l|. The equilibrium coexistence for 
given overall composition ($^,^^) minimises the total 
free energy F subject to constraints relating the bilayer 
phases’ area fractions On and compositions to the 
overall compositions: 

phases phases phases 

0n = l,Yl , E • (») 

F is proportional to the height of the tangent plane at 
the given (^^,^^). Coexistences not fully minimising F 
are metastable. 

Partial phase diagrams for two of the main parame¬ 
ter points to be used in this work (except Eigs. |H1 fTTl 
and Section IVDj) appear in Eig. [2)3, c. Spinodal lines en¬ 
close regions of local stability. An overall composition 
(<|)t^<|)b^ within a tie-line (two-phase) or triangle (three- 
phase) can split into coexisting bilayer phases defined by 
the endpoints. The bilayer phases are, qualitatively, R 
(both leaflets dominated by the same species, approxi¬ 
mate compositional symmetry between leaflets) or AR 
(opposite species, strong compositional asymmetry). R 
phases thus contain a majority of either SS or UU sites, 
and AR mostly SU or US. The precise composition of a 
given phase depends quantitatively on which phase coex¬ 
istence it is part of. The separate leaffet compositions in 
a bilaw phase are given by the projections onto the 
axes [2I, y, [l0| . A state of phase coexistence is denoted 
R-R (two compositionally symmetric bilayer phases), R- 
R-AR (two symmetric phases and a strongly asymmetric 
one it AR-AR [H, etc. Overlapping tie-lines or tri¬ 
angles on the phase diagram imply that multole phase 
coexistences are possible for a given (^^, ^^) [l|. Eig. [2^ 
establishes shorthand diagrams for the equilibrium and 
metastable coexistences appearing in Eig. [2)3,0. 

In 0 it is claimed that, for overall symmetric leaflets 
(^b = “perfect antiregistration” is possible only if 
^b_^t_Q 5 (cf. Eq. [13]). It is true that AR-AR coexis¬ 
tence (strong local asymmetry throughout the bilayer) is 
possible only in a limited region centred on = <I>^ = 0.5 
(see Appendix EJ. Similarly, conventional R-R coexis¬ 
tence is available for <|>^ = <!>*, but may not be for 
(in [ 2 I, this lead to R-R-AR coexistence, corresponding 
to the black triangles on Eig. [2)3,0). However, even for 
overall leaffet compositions that prohibit one or both of 
“perfect antiregistration” (AR-AR) and conventional R- 
R coexistence, Eig. [ 2 ] shows that there can still be mul¬ 
tiple available phase eoexistenees. The available states 
contain greater or lesser amounts of R versus AR phases 
and hence degrees of symmetry than the hypothetical 
case of fully independent leaffets. 

Eor instance, AR-AR-R coexistence can be considered 
the “minimally registered” state where overall composi¬ 
tion prohibits AR-AR coexistence; all bilayer regions are 
antiregistered except for one residual R phase. AR-AR-R 
may compete with the more registered R-R (for <I>^ = ^^) 
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FIG. 2. (Color online) (a) Shorthand diagrams for the metastable (AR-AR, AR-AR-R) and equilibrium (R-R, R-R-AR) 
coexistences studied in this work. (b,c) Grey contour plots of free-energy density /(0^, 0^) with partial phase diagrams overlaid. 
Ao — 2a^ Hi — 3a~^kBT, V = 0.6/cbF, with J and B as indicated. Two-phase central R-R and AR-AR tie-lines (dashed) and 
three-phase triangles (solid or dotted) are shown. Metastable coexistences are red. A comprehensive version of (b) and detailed 
discussion appears in Appendix [A] 


or R-R-AR [H (for <|>^ 7 ^ <l>^) states. This behaviour is 
shown in the simulation results presented below. 

As explained in an upcoming Comment [lH, the AR- 
AR-R state is relevant in evaluating the claim of [HI 
that common observations of R-R can be explained via 
line energies alone, without a direct inter-leafiet coupling 
that would lead to an explanation from bulk free energies. 
In the Comment [28| we show that the elastic theory in 
[HI implies that R-R may be lower in total line energy 
than AR-AR-R for some compositions and parameters, 
but that for other compositions and parameters R-R is 
not equilibrium on the basis of line energies alone. 

Fig. [ 213 , c are qualitatively similar; in both (as in the 
schematic landscape Fig.[T]) R minima are lower than AR 
so that equilibrium coexistence is R-R or R-R-AR (or 
R-AR, see Appendix [3, depending on overall composi¬ 
tion. A larger direct coupling B (Fig. [JJd) further raises 
the free energy of AR minima due to the energy cost of 
SU and US sites, hence the phase coexistence regions 
are quantitatively different between Fig. [ 2 l 3 ,c. Other pa¬ 
rameter choices can yield AR minima lower than the R, 
so that AR-AR and AR-AR-R become equilibrium, not 
metastable, states. 

For simplicity, Fig.[2]shows only the tie-lines and trian¬ 
gles studied in the present simulations. A comprehensive 
version of Fig.[ 2 l 3 , with further discussion of bilayer phase 
diagrams in the context of existing work and experimen¬ 
tal data, appears in Appendix [Al 


C. Registration is nonconserved 

Since we assume no flip-flop or solvent exchange, the 
overall leaflet compositions are conserved [9|. How¬ 
ever, depending on the phase coexistence the bilayer 


chooses, the overall number of pairwise registered lipids 
varies; AR {SU or US) and R {SS or UU) sites inter¬ 
convert via 

su^us^ss^uu. 


The nonconserved overall occupancies Afa of each site 
type a G {SS, UU, SU, US} are constrained by the con¬ 
served overall leaflet compositions 


Msu - J^us = , 

(9) 

^fss + ^fsu = ^W , 

( 10 ) 

J^uu + Mus = (1 - . 

( 11 ) 

This leaves a free variable 


. Vss + Afuu 

M 

( 12 ) 


the degree of molecule-level transbilayer symmetry - the 
proportion of lipids, over the whole bilayer, that micro¬ 
scopically appose one of the same species. It can vary in 
the range 

- 1| < A < 1 - - $'"1 . (13) 

Hence, as well as domain coarsening of a given coexis¬ 
tence of phases (during which A is roughly constant), 
the bilayer may switch between competing phase coex¬ 
istences (Section [HB[) . leading to a significant change in 
A while still conserving the overall leaflet compositions 
^t(b) Therefore, A enables us to monitor such transi¬ 
tions. Its value Ao at t = 0 corresponds to the high- 
temperature limit: 

Ao = + • 


(14) 
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Ao can be interpreted as the value that would be taken if 
the leaflets were completely independent of one another, 
coupled neither by a direct interaction or their combined 
hydrophobic thickness. 

Experimental measurement of A could be achieved us¬ 
ing FRET (Forster resonance energy transfer), with the 
donors and acceptors distributed exclusively in opposite 
leaflets. 


D. Qualitative mapping onto conventional ternary 
phase diagrams 

To properly account for transbilayer symme¬ 
try/asymmetry in describing phase behaviour, our 
model uses a separate composition axis for each leaflet 
H i We now discuss how conventional bilayer 

phase diagrams map onto this description. 

Fig. [3^ shows a typical lipid phase diagram of a ternary 
bilayer, which is taken to describe the fully registered case 
in which exact compositional symmetry holds through¬ 
out the bilayer. The binary cholesterol-free edge (sat.- 
unsat.) of Fig. [3^ corresponds to the 0^ = 0^ diagonal of 
Fig. [3b. The segment of that edge linking Ld and gel cor¬ 
responds to the R-R tie-line of Fig.|3b (squares). Because 
the leaflet compositions of R phases involved in R-R- 
AR coexistence will, in general, differ quantitatively from 
those involved in R-R coexistence, the mapping to a con¬ 
ventional ternary phase diagram is non-trivial. R phases 
in the R-R-AR region (circles and triangles) of Fig. |3b 
do not he precisely on the (j)^ = (j)^ diagonal, so only 
have approximately symmetric local leaflet compositions. 
Thus, for the corresponding symbols on Fig. |3b the top 
and bottom leaflet compositions are slightly separated. 
The leaflet compositions for the strongly-asymmetric AR 
phase (open circle) lie towards opposite ends of the T^^-gel 
binodal. 

Similarly, a binary model can represent a ternary bi¬ 
layer in, e.g., an Lg-Ld reg ion of Fig. [3^. This idea is 
exploited in simulation |29[ 1^ and theory i i i, [ii , 
and a careful derivation by Garbes Putzel and Schick @ 
shows that it entails assuming the key compositional or¬ 
der parameter is the relative abundance of saturated ver¬ 
sus unsaturated lipids. Two-component Ising universal¬ 
ity in Lg-Ld experiments [3l| supports such a mapping. 
Fig. Eb illustrates this mapping for our bilayer phase di¬ 
agram. The R-R tie-line of Fig.|3t (squares) corresponds 
to an Lg-Ld tie-line of Fig.|3b, but the phases involved in 
the R-R-AR coexistence region need not map to leaflet 
compositions on precisely the same Lg-Ld tie-line. Thus, 
the correspondence of a binary model to a pseudo-binary 
region of a real ternary system is approximate but the 
important qualitative features should be preserved. 



FIG. 3. (Color online) (a) Schematic conventional ternary 
phase diagram, (b) Leaflet compositions for coexistences in 
a binary (lipid-lipid) bilayer, using R-R and R-R-AR coex¬ 
istences as examples. Pink symbols in (b) correspond to 
pink symbols in (a), and letters “t” and “b” in (a) label 
whether the symbol refers to top or bottom leaflets, or both, 
(c) For pseudo-binary Lg-Ld coexistence, the relationship to 
the ternary phase diagram is less precise. The corresponding 
blue symbols in (a) have halos indicating the uncertainty in 
ternary phase space. 


III. SIMULATION PROTOCOL AND 
PARAMETERS 

A. Protocol 

We simulate the model with a Kinetic Monte Garlo al¬ 
gorithm resembling Kawasaki (spin-exchange) dynamics 
1^ . The moves are: lateral exchanges of lipids within a 
single leaflet; coordinated lateral exchanges, where the 
top and bottom lipids of a lattice site move as one; 
and hydrophobic length changes of single lipids, pro¬ 
posed uniformly between ±0.5 a. The total lattice size 
is = 200^ = M. Time t is measured Monte Garlo 
Steps (MGS), comprising 2M lateral exchanges, 2M com¬ 
bined lateral exchanges, and 2Mnr length changes. We 
set fir = 10, so that lipids relax their lengths faster than 
they laterally diffuse Q. The bilayer is initialised with 
the lipids in each leaflet randomly distributed and the 
lengths initially relaxed ^ 1000 times per lipid. The sim¬ 
ulation models an instantaneous quench from high tem¬ 
perature at t = 0. 

With idealised dynamics it is difficult to assign a cor¬ 
responding physical time unit to the MGS. One could 
attempt to use characteristic self or collective diffusion 
times of the simulation, but a more straightforward com¬ 
parison to experiment or molecular simulations could be 
made on the basis of domain size after a queneh^ rather 
than time after a quench. We can note that of is the lat¬ 
eral area of one lipid. (Thus, for example, the crossover 
to registered domains occuring in Fig. [6] is complete by 
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the time the domain lengthscale reaches ^ 100 a.) In 
we calculate nucleation energetics to estimate criti¬ 
cal domain sizes for registration, presumably related to 
the lengthscale at which this crossover takes place. 

Similarly, hydrodynamics are not included, but will 
dominate the kinetics as domains exceed a certain length 
scale Q, estimated by Fan et al. [35| to be ^ 10“^ m 
in lipid bilayers. However, hydrodynamics cannot affect 
the free-energy landscape, which determines the compet¬ 
ing metastable and equilibrium states. Hence we expect 
the competition of registered and antiregistered phases 
studied here to remain robust. It is of great interest to 
consider marrying hydrodynamics to a coarse-grained bi¬ 
layer free energy that accounts for both leaflets and their 
couplings, such as that derived from our model [l|. 


B. Registered and antiregistered structure factors 


We also monitor domain growth using separate struc¬ 
ture factors for R and AR lattice sites (where 

m = R, AR). Their first moments [36| 


'EqQSmiq^t) 


(15) 


lead to characteristic lengthscales for domains of micro¬ 
scopically registered versus antiregistered lipids 


Sm{t) 


27r 

qm{t) 


( 16 ) 


In the results presented here Sm{t) is averaged over 
twenty independent trajectories. S^{q^t) can be related 
to more conventional structure factors suitable for exper¬ 
imental measurement (see Appendix [B]) . 


C. Parameters 

Detailed discussion of parameterisation is given in [l| . 
The lattice spacing is a ^ 0.8 nm. We use k = 3a~‘^kBT 
{ka ~ GOteTum”^, in the range for lipid bilayers at 
300K a fiducial estimate of the indirect cou¬ 
pling parameter is J ^ 2a Q. B leads to an effec¬ 

tive value of 7 which has been wid^ estimated in the lit¬ 
erature as 7 ^ 0.01 — I kBTnm~‘^ [El, [I, QIH, 113, [3 • In 
most simulations here we set the length/structural mis¬ 
match parameter Aq = 2a ^ 1.6nm. This relatively 
large value strengthens the indirect and direct coupling 
effects, making the kinetic regimes in the simulations 
clearer to interpret. Simulations with a smaller value 
Ao ~ I a yield similar kinetics fSection IV Pp . i.e., the 
balance of J and B determines what kinetics occurs. 

In the theoretical plots (Figs. [2] and 0]) we use V = 
0.6 /cbT to exceed the threshold Vq = 0.5 /cbT required, 
in mean-field, for phase separation in the absence of any 
other couplings. However, the equivalent threshold in 
simulation (where fluctuations are included) is = 



coupling / k^T 


Aa; arb. units 



FIG. 4. (Color online) Stability diagram for a bilayer of 
(0.5, 0.5) overall composition, with Ao = 2 a, a = Sa^^/csT, 
V — 0.6/cbT, calculated from the mean-field free energy [I[. 
The equilibrium is R-R coexistence (R free-energy minima are 
lower than AR) except for a tiny region BcB jk^T < 0.005 
which we do not study here. Below the R instability line, the 
initial homogeneous state is not unstable to the R mode al¬ 
though R-R separation is still equilibrium [I|. Above the AR 
instability line, the AR minima do not exist (Appendix 
and there is no instability to the AR mode. Colours show the 
relative growth rates of R and AR modes. 


0.88 k^T. We employ V = 0.9 k^T for the simulations in 
an attempt to ensure a similar qualitative regime as the 
theoretical predictions. We have also tested the theory’s 
sensitivity to V by replotting Fig. 0] with V = 0.9 ksT 
(Fig. [T7| in Appendix El, finding only a small difference 
in the predicted R/AR competition. This suggests the 
precise value of V is not of crucial importance to our 
present study, so long as R > Vb (mean-field) or R > 
(simulation). 


D. Flip-flop 


The present modelling does not include flip-flop. 
Wholesale redistribution of lipids between leaflets may 
influence phase behaviour [3, lio| . In 0 , bilayers of asym¬ 
metric overall leaflet compositions ^ appeared to 
relax over hours to overall symmetry through 

flip-flop. Such progression toward overall symmetry is 
perhaps the most intuitive effect of flip-flop, but in j^ . 
in contrast, a slow decay towards an overall asymmetric 
state was inferred. While it seems that experiments on 
phase separation can be safely performed before flip-flop 
becomes important Sin, we are currently extending 
our modelling to include the effects of flip-flop. 
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IV. KINETIC CONSIDERATIONS 
A. Initial instability and nucleation kinetics 

Consider ^^ = ^^ = 0.5, which we label (0.5, 0.5). After 
quenching to an unstable region, the local composition is 
initially uniform at (0.5, 0.5). Linear stability analysis 
then yields growth rates of symmetric (R) or asymmet¬ 
ric (AR) modes corresponding to demixing perturbations 
along the R-R or AR-AR central tie lines of Fig. (Jb [l|, 
illustrated by the curved arrows on Fig. [TJ The analysis 
involves a Ginzburg-Landau free energy 

FG.L = y'dV((//a2) + /g,,d) , (17) 

where/grad = \J{Vdf . The bulk 

term / drives phase separation, while /grad penalises re¬ 
sultant gradients in composition and bilayer thickness. 
The balance determines which mode is faster and thus 
dominates initial phase separation (Fi^|4]). Further de¬ 
tails can be found in Appendix O and ^ . 

For both R and AR modes, the Ising-like term (length- 
independent interactions, V) promotes instability via / 
and penalises composition gradients via /grad- The 
drophobic mismatch penalty J plays a subtler role [i[. 
In /grad it penalises thickness gradients (only incurred 
by the R mode), without necessarily boosting R instabil¬ 
ity through the bulk term /. Hence, the AR mode can 
be faster although R-R coexistence is equilibrium. This 
is an explicit derivation of the physics underlying Ost- 
wald’s heuristic rule of stages (in this particular system), 
i.e., that a system will access metastable states on its 
way to the lowest free energy state. Away from (0.5, 0.5) 
overall composition, competing instabilities will be more 
complex than the modes studied in Fig. 01 but the com¬ 
petition between local symmetry and local asymmetry 
remains similar. Parameter ranges for typical phospho¬ 
lipids imply that these competing modes will be of com¬ 
parable magnitude [l| so that, e.g., increasing the tail 
length mismatch could tip the balance to make the AR 
mode dominant. 

Nucleation is required to reach equilibrium if initial 
demixing leads to a metastable state, although a bilayer 
of (0.5, 0.5) composition is conventionally thought to sep¬ 
arate directly into registered phases by spinodal decom¬ 
position. In the cases studied here, the metastable coex¬ 
istence is AR-AR or AR-AR-R, and equilibrium is R-R or 
R-R-AR, thus requiring nucleation of registered phases. 
In the nucleation process an area-dependent energetic 
benefit of converting AR to R phases competes with a 
cost from thickness mismatch at the R domain bound¬ 
ary; this is studied in H. The hydrophobic mismatch 
penalty J penalises thickness mismatch at the bound¬ 
ary of a nucleating registered domain. Thus, as well as 
kinetically favouring AR-AR separation immediately af¬ 
ter the quench, hydrophobic mismatch inhibits the nu¬ 
cleation required to equilibrate from a metastable state. 



HUU RSS 


FIG. 5. (Color online) Schematic free-energy landscape of 
Fig. m viewed in the (—1,1) direction (arrow, inset), such 
that AR minima appear one in front of the other. We super¬ 
impose tie-lines and triangles corresponding to the multiple 
coexistences available for (0.5, 0.5) overall composition: R-R 
equilibrium (black); two metastable AR-AR-R triangles (lines 
from this angle, red); metastable AR-AR (dot from this an¬ 
gle, red). One of the AR-AR-R triangles is illustrated on the 
inset. 


The direct coupling B counters both these effects. It 
promotes initial instability to the R mode, and raises the 
AR free-energy minima relative to R (Fig. [2]3 versus c), 
so increasing the area-dependent benefit of nucleating the 
equilibrium R phases. 

Thus, depending on the competing inter-leaflet cou¬ 
plings, three classes of kinetics are possible H; direct 
separation into equilibrium phases; nucleation out of a 
metastable state; and trapping in a metastable state. 


B. Multiple metastable states 

A (0.5,0.5) bilayer can he within an R-R tie-line, two 
metastable AR-AR-R triangles and a metastable AR-AR 
tie-line, and thus access multiple metastable phase co¬ 
existences. As illustrated for the schematic free energy 
landscape in Fig. [5l the total bulk free energy at AR- 
AR-R coexistence - given by the height of the tangent 
plane at (0.5, 0.5) - is lower than for the AR-AR cen¬ 
tral tie-line. Thus, AR-AR coexistence can convert to 
AR-AR-R for a slight drop in bulk free energy (or, of 
course, reach equilibrium R-R to fully minimise free en¬ 
ergy). The AR-AR-R planes are degenerate for (0.5, 0.5), 
so a metastably trapped bilayer could locally fluctuate 
between SS or UU R phases in different regions (this 
seems to occur in Section IV PI where weaker inter-leaflet 
couplings make fluctuations more apparent). In contrast, 
an off-equimolar overall composition on the (j)^ = (j)^ diag¬ 
onal (e.g., (0.3,0.3)) can access a single AR-AR-R state, 
though equilibrium is still R-R. Appendix contains a 
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detailed diagram of all coexistence regions from the phase 
diagram in Fig.lJJ), and discussion of their interpretation 
in relation to experiments. 


V. SIMULATED PHASE-TRANSITION 
KINETICS 

We now use simulation to investigate the kinet¬ 
ics of competing phase coexistences. We vary the 
overall leaflet compositions to study different 

metastable/equilibrium coexistences, and vary the indi¬ 
rect (hydrophobic mismatch) coupling J and direct cou¬ 
pling B to explore different kinetic regimes which we ex¬ 
pect to lie within the parameter ranges of typical phos¬ 
pholipids. Increasing J could be physically realised by 
increasing the tail length of the saturated lipid species to 
increase hydrophobic mismatch [a. Increasing B might 
be achieved by increased difference in tail unsaturation 
between species, increasing their structural mismatch. 
Unless otherwise noted, the other simulation parameters 
are Aq = 2 a, tv = 3a~‘^kBT, V = 0.9/cbT. As well 
as visual inspection, phase-separation kinetics are moni¬ 
tored by the degree of microscopic transbilayer symme¬ 
try A (which could be experimentally accessed, e.g. . by 
FRET). Visualisation was performed with OVITO (i^ . 
Videos corresponding to Figs. IbHTOl are available online 
0 . 


A. ($*, $’^) = (0.5,0.5): R-R versus AR-AR 

For (0.5,0.5) overall composition, metastable AR-AR 
coexistence (local asymmetry throughout) competes with 
R-R (local symmetry throughout). For J = 4a“^/cBT, 
B = OAS a~^kBT (Fig. [6j), we observe two-step kinet¬ 
ics. There is an immediate significant drop in A because 
J > B favours pairwise antiregistration at the single¬ 
lipid level in the laterally homogeneous initial state [l| 
(in agreement with the experiments of HH). Metastable 
AR domains form and coarsen (Fig. [Ob). The bilayer re¬ 
duces the hydrophobic cost of remaining R sites through 
annihilation into further AR sites, and coalescence into 
R domains. These coalesced groups become nuclei of SS 
and UU phases [23|. From t ^ 10^ — 10^, these nuclei 
grow at the expense of the AR phases, increasing A and 
converting the bilayer from metastable AR-AR to equi¬ 
librium R-R coexistence for a payoff in bulk free energy. 
The stochastic nature of the nucleation transition is ap¬ 
parent by the spread in A between independent simula¬ 
tion trajectories during this period. Thereafter, coarsen¬ 
ing of the equilibrium phases continues and A is roughly 
constant. 

Hence, the transition between the competing phase co¬ 
existences is signified by nonmonotonic evolution of A. 
The bilayer first becomes more asymmetric than the un¬ 
correlated case (A = Aq) as it favours AR-AR coexis¬ 
tence, then more symmetric as it switches to equilib- 


metastable metastable 




FIG. 6. (a) Grey: degree of pairwise lipid symmetry A 

through time (twenty independent runs) for a simulated bi¬ 
layer of (0.5, 0.5) overall composition with J = Aa~^kBT, 
B = 0.48a“^A;BT, Ao = 2 a, /€ = Sa^^/csT, V — 0.9/cbT. 
The initial (uncorrelated leaflets) value is Ao = 0.5. Black 
(secondary axis): characteristic lengthscale of R and AR do¬ 
mains (averaged over the twenty runs), (b) Trajectory snap¬ 
shots showing the transition from AR-AR to R-R coexistence, 
causing nonmonotonic variation of A. For the second row of 
snapshots, diagrams of the competing coexistences are shown, 
with overall composition marked by a dot. Videos correspond¬ 
ing to Figs. [BHTOl are available online 


rium R-R. The two-step kinetics in Fig.[6]is in agreement 
with Fig. m J and B for this parameter point are well 
within the region where the mean-field theory predicts 
the AR mode to dominate so that AR-AR occurs first, 
as in Fig. [6l 

The lengthscales 5 r and sar show initially faster 
growth of AR domains. During conversion of AR to R 
domains, there is an increase in the growth exponent of 
registered domains to q^r 0.5 (where sr ^ t^^), while 
AR domain growth slows. This exponent is expected for 
the nonconserved Ising (Model A) universality class H, 
and we attribute it to nonconserved growth of R nuclei 
via annihilation of AR domains. As A begins to satu¬ 
rate, further conversion of AR slows. Thus, SS and UU 
now play the role of “up” and “down” spins in the con¬ 
served Ising model. R domain growth slows to q^r 0.28, 
a typical apparent exponent for diffusive coarsening in 
the conserved Ising (Model B) class [36|, which would 
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FIG. 7. J = Aa-‘^k^T, B = 0.24 (0.5, 0.5) overall 

composition. As Fig. [6l with weaker direct coupling B. Nu- 
cleation of R domains is suppressed and the bilayer is trapped 
in AR-AR coexistence. 


FIG. 9. J = 4a“^/cBT, B = 0.48a“^/cBT, (0.3, 0.3) overall 
composition. As Fig.[6l for (0.3, 0.3) overall composition. The 
bilayer moves from AR-AR-R to R-R coexistence. 



FIG. 8. J = OAa-^kBT, B = 0.48 (0.5, 0.5) overall 

composition. As Fig.[6l with weaker indirect coupling J. The 
R mode is fastest, and the bilayer immediately accesses R-R 
coexistence. 


tend to a = -^ in the asymptotic time limit. We employ 
a ~ 0.28 as a heuristic test for tending towards the con¬ 
served growth scaling regime. As domains grow larger 
these exponents will be altered by hydrodynamics (Sec¬ 
tion iTTm . which is not included in our model. 

With weaker direct coupling B = 0.24a“^/cBT (Fig.[7|) 


the bilayer gets trapped in AR-AR coexistence and all 
R nuclei are destroyed, leaving the metastable phases to 
coarsen indefinitely. AR domains coarsen with an expo¬ 
nent q^ar ~ 0.28 reflecting conserved growth, SU and 
US now playing the role of up and down spins in the 
conserved Ising model. This metastably trapped AR-AR 
state resembles that found in molecular simulations with 
strong ^drophobic length mismatch [mil, although 
in I 12 L 1 131 it could be stabilised due to finite simulation 
size 1^ . 

For B = 0ASa~‘^kBT and weaker hydrophobic mis¬ 
match parameter J = 0.4a“^/cBT (Fig. [8]) the bilayer 
proceeds straight to R-R coexistence. The proportion of 
registered lipids A increases immediately, in contrast to 
the two-step behaviour seen in Fig.O R domains reach a 
growth exponent q^r ^ 0.28 without a regime of noncon- 
served (q^r ~ 0.5) growth. This behaviour is consistent 
with the prediction of Fig. 0] that for these parameters 
the R mode is faster during initial phase separation. 


B. (0.3, 0.3): R-R versus AR-AR-R 

A common experimental case is leaflets containing the 
same, but nonequimolar mixture. If the composition is 
within the R-R central tie-line, equilibrium coexistence 
is R-R. We use = 0.3, entailing excess UU sites 

relative to the (0.5, 0.5) case. According to Figs. [2]3,c, it 
is also within a metastable AR-AR-R coexistence region. 

Fig.[9]shows the associated two-step kinetics. It is anal¬ 
ogous to that for (0.5, 0.5) composition in Fig. [6l except 
that now the metastable state comprises two AR phases 
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metastable 


metastable 



FIG. 10. J = 4a-^kBT, B = OASa-^ksT, (0.5, 0.7) overall 
composition. As Fig. El for (0.5, 0.7) overall composition. The 
bilayer moves from AR-AR-R to R-R-AR coexistence. 



FIG. 11. Trajectories using Aq = la, which weakens both 
indirect and direct couplings [I[. The overall composition is 
(0.5,0.5). B = 0.3 a~^kBT, and J is increased, (a) J = 
O.la^^/csT. Immediate separation into R-R coexistence (cf. 
Fig. [8|). (b) J = la~^kBT. Competing AR-AR and R-R 

coexistences, slowly converting to R-R domination (cf. Fig. [6]). 
(c) J — Aa~^kBT. Metastably trapped in AR-AR (cf. Fig.[7|). 


plus a UU R phase. Again, A initially drops from the un¬ 
correlated case Ao, this time indicating a preference for 
the AR-AR-R state, then increases, as R-R is reached. 
Unlike in Fig. [6l there is no clear ckr ^ 0.5 growth dur¬ 
ing AR to R conversion. This probably arises because the 
metastable state now contains AR-AR-R, not just AR- 
AR, so the transition to equilibrium coexistence involves 
not only AR to R conversion but also continued coars¬ 
ening of the existing UU domains. This is also shown in 
the behaviour of A which, in comparison to Fig. [6l does 
not attain as low a minimum, and so grows less dra¬ 
matically. With smaller B (Appendix [D]), nucleation of 
the registered SS phase is suppressed, leaving metastable 
AR-AR-R coexistence. 


These simulations thus demonstrate an important con¬ 
sequence of the discussion in Section ing even where per¬ 
fect antiregistration (AR-AR) is not possible, “minimally 
registered” AR-AR-R is possible. This is in contradiction 
to , which uses the fact that AR-AR is only possible 
in a narrow window of compositions (Fig. IT^) to argue 
that R-R is generally stabilised by line energies alone, so 
that explaining R-R would not require direct inter-leafiet 
coupling (a bulk free energy effect). In an upcoming Com¬ 
ment we show that R-R is stabilised by line energies 
only for certain compositions and parameters, because 
one must compare it to the AR-AR-R state if AR-AR is 
not available [ 43 . 


C. (0.5, 0.7): R-R-AR versus AR-AR-R 

We now consider asymmetric overall compositions, as 
in experiments on bilayers specially prepared with differ¬ 
ing leaflets [2|. Overall composition (0.5, 0.7), according 
to Figs. [2)3,0, gives R-R-AR coexistence at equilibrium, 
as was observed in [2[. Just as (0.3, 0.3) prevented per¬ 
fectly asymmetric AR-AR, here (0.5, 0.7) prevents per¬ 
fectly symmetric R-R. 

This R-R-AR state competes with AR-AR-R coexis¬ 
tence. In Fig. uni two-step kinetics can be seen. We 
effectively have coarsening of a 2D ternary mixture {SU, 
US and 7/7/-dominated phases), as recently studied in 
[ 4 ^ . Accordingly, we expect the same scaling as for a 
binary mixture and, indeed, both measures of domain 
size exhibit ckr ~ ckar ~ 0.28 growth. With smaller 
B (Appendix iD]) . the bilayer is metastably trapped in 
AR-AR-R. If the overall leaflet compositions are strongly 
asymmetric, the equilibrium state may be R-AR instead 
of R-R-AR. This is discussed in Appendix in relation 
to some available experimental findings [I, |l6|, |4l| . 


D. Weaker mismatch 

So far we have used a mismatch parameter Aq = 2 a ~ 
1.6 nm to strengthen the indirect and direct couplings 
and more easily distinguish the kinetic regimes. With 
a smaller value Aq = la, both the indirect and direct 
couplings become weaker, but estimated phospholipid 
parameters still imply competing R and AR modes of 
































































11 



FIG. 12. (Color online) (a) Snapshots from the trajectory from Fig. [9] rendered in pseudo-3D. Colours showing total bilayer 
thickness as would be measured via AFM (analogous to total bilayer fluorescence in standard fluorescence microscopy) are 
compared with greyscale showing leaflet registration and antiregistration (as in Fig. [9]). (b) Snapshots from a trajectory 

with weaker direct inter-leaflet coupling B (see Appendix [DJ. Nucleation of SS is suppressed, leaving metastable three-phase 
coexistence which could be mistakenly interpreted as two-phase equilibrium in the thickness representation. 


comparable magnitude [l|. Fig. [TT] shows the same phe¬ 
nomenology with Ao = 1 a as in the previous simulations. 
As expected, increasing indirect coupling J in Fig.fTTlhas 
a similar effect to decreasing B: for small J the R mode is 
fastest; for intermediate J, AR and R compete; for larger 
J, metastable trapping in AR-AR coexistence is evident 
as the energy cost for thickness mismatch prevents nu¬ 
cleation of R phases. In Fig. [TTh . an AR-AR trapped 
state appears to locally fluctuate into AR-AR-R, creat¬ 
ing transient regions of R phase. This is consistent with 
the discussion in Section ITVBl 


VI. IMAGING THAT DOES NOT RESOLVE 
DISTINCT LEAFLETS 

Phase-separating bilayers are often imaged with the 
fluorescent tag distributed throughout both leaflets. In 
this method, one does not measure the separate leaflets, 
but rather the sum of the two leaflets’ compositions. 
This is analogous to measuring bilayer thickness via AFM 
(atomic force microscopy), by which the sum of the two 
leaflet thickness is measured - coexisting AR phases {SU 
and US) would be of the same thickness and thus not 
distinguishable from one another. We imitate this by 
imaging the simulated bilayer in a total thickness repre¬ 
sentation (Fig. [12]), no longer labelling SS^ SU etc. The 
bilayer is rendered in pseudo-3D by assigning total thick¬ 
ness to 2 ;. (The reference do = Iso S- ^uo is irrelevant to 
the behaviour of the model, but for Fig. [12] we have set 
do = 5 a ^ 4 nm.) 

Fig. [12^ shows snapshots corresponding to Fig. [9] Dur¬ 
ing the transition from metastable to equilibrium co¬ 
existence, three distinct thicknesses, arising from four 
distinct phases, appear. The bilayer initially prefers 


metastable AR-AR-R coexistence, where the R phase is 
UU. The two AR phases have identical, intermediate 
thickness, UU is thinner and small nuclei of SS are thick¬ 
est. The SS and UU phases then grow at the expense of 
AR {SU and US) phases, leaving the two thicknesses of 
the SS and UU phases. The appearance of three distinct 
thicknesses and nonconservation through time of their 
respective area fractions would thus signify a transition 
between competing coexistences in molecular simulation 
or experiment. AFM would be well-suited, so that these 
kinetics could be monitored while domains are too small 
to resolve optically. 

In Fig. [T2]3 the bilayer remains trapped in metastable 
AR-AR-R coexistence. In the absence of an SS phase for 
comparison (or a priori knowledge of its expected thick¬ 
ness) it is in principle possible to misinterpret the two AR 
phases as a single, registered phase. Hence, long-lived or 
steady state (i.e., trapped) metastable three-phase coex¬ 
istence could masquerade as two-phase (R-R) equilibrium 
coexistence. The apparent area fractions of thicker and 
thinner phases differ from the true R-R equilibrium state. 

However, studies in which leaflets are independently 
fluorescently labelled suggest the expected equilibrium 
of R-R-AR Q or R-R coexistence [ij is reached, imply¬ 
ing that such metastable trapping is avoided. Therefore, 
it is likely that experiments with similar model phospho¬ 
lipids also reach equilibrium phase coexistence, and the 
misinterpretation shown in Fig. [T^ . requiring long-lived 
metastability that is apparently equilibrium, does not 
arise. However, experiments limited by optical resolution 
would not probe the kinetics of a metastable to equi¬ 
librium transition on time/lengthscales before domains 
reach resolvable size. Large hydrophobic mismatch would 
increase the possibility of long-lived metastable AR-AR 
as seen in simulation [ill, [ 2 ^, while Lin et ah [io| found 
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slow creation of AR domains at the expense of R, in a 
liquid-gel bilayer (perhaps influenced by the solid sub¬ 
strate). 

Labelling the leaflets separately provides extra insight 
by revealing whether the pattern of enrichment and de¬ 
pletion of dye in one leaflet is colocalised with that in 
the other. However, such colocalisation need not neces¬ 
sarily imply R-R coexistence (local compositional sym¬ 
metry everywhere), but can also indicate R-AR (see Ap¬ 
pendix and Q). 


VII. DISCUSSION 

Evidence in the literature implies that the leaflets of 
a phase-separating bilayer are inherently coupled, both 
by a direct coupling {B) that encourages compositionally 
symmetric bilayer phases (registered, R), and an indirect 
coupling (J) from hydrophobic mismatch that favours 
strong compositional asymmetry (antiregistered phases, 
AR) [i3, [3 [3, 134^ . Eoj . Here we have studied the 

resulting kinetics, by simulating the lattice model un¬ 
derlying the theory in [l|. We have found signatures by 
which the kinetics could be detected in molecular simu¬ 
lation or experiment, to provide insight into fundamental 
interactions. 

Governed by the inter-leaflet couplings and the over¬ 
all leaflet compositions, multiple phase coexistences in¬ 
volving different degrees of transbilayer symmetry can 
compete. An important result of the present work is 
that this applies over wide regions of the phase dia¬ 
gram. For instance, AR-AR coexistence (i.e., strong lo¬ 
cal asymmetry everywhere throughout the bilayer, which 
has been called “perfect antiregistration” 0 ), may only 
be possible in a small window of overall compositions 
(see Appendix El). However, outside this window, “min¬ 
imally registered” AR-AR-R coexistence can still occur, 
as well as the typically-observed R-R Q or R-R-AR Q 
states dil. Hence, for a variety of overall compositions, a 
phase separating bilayer can choose between states that 
are either more or less symmetric than the hypothetical 
“uncorrelated” case of completely independent leaflets 
(Eq. [14]). 

The direct inter-leaflet coupling interaction B typi¬ 
cally renders the more-asymmetric (AR-AR or AR-AR- 
R) states metastahle^ due to their higher bulk free energy, 
providing an appealing explanation for observations of 
domain registration. Conversely, observed domain an¬ 
tiregistration [mil may be understood as metastably 
trapped due to strong hydrophobic mismatch, or sta¬ 
bilised by small domain size • If was argued recently 
that domain registration observations are explained in 
general by line energies alone, without direct inter-leaflet 
coupling [TqI - our reasons for disputing this claim are 
detailed in an upcoming Comment [28|. 

Physical parameter estimates [l| imply that competing 
instabilities to form registered versus antiregistered do¬ 
mains can be of comparable strength, so that changes to 


molecule tail length or structure mismatch - effectively 
tuning J and B - may determine which occurs first af¬ 
ter a quench. Subsequent to reaching a metastable state, 
the bilayer must nucleate to reach equilibrium We 
also find a preference for pairwise lipid asymmetry in the 
initial homogeneous state, in line with Hence, in¬ 

cluding both direct inter-leaflet coupling and hydropho¬ 
bic mismatch effects helps unify observations of trans¬ 
bilayer symmetry and asymmetry [H MM, 0 

in mixed bilayers. The competing inter-leaflet couplings 
also lead to a “critical radius” for domain registration, 
below which AR domains are stable [3; fhis could be 
important in cellular rafts or clusters, which are thought 
to exist in a size range 10 — 100nm [49| [5^ . 

The findings here can be probed by studying the kinet¬ 
ics of transbilayer symmetry upon quenching from high 
temperature. Conducting such studies either in molec¬ 
ular simulations or experiment can, via our theory, di¬ 
rectly reveal the relative importance of the direct versus 
indirect inter-leaflet couplings [l|. Therefore, although 
the simulation used here is idealised fSection IHI Ap . it 
has the benefit of directly linking observed kinetics to a 
free-energy landscape derived from the same microscopic 
model; this is turn provides a basis for more realistic sim¬ 
ulations and experiments to test the essential features 
captured in the model. An interesting experiment would 
be a FRET study in which donors and acceptors are dis¬ 
tributed in opposite leaflets, thus accessing the degree 
of microscopic transbilayer symmetry A. A nonmono¬ 
tonic evolution of the FRET signal after a rapid quench 
to a phase-separating region would be a clear signature 
of competing asymmetric and symmetric modes of phase 
separation. 


Appendix A: Bilayer phase diagram details 

Fig. [13] provides detail on the phase diagram from 
Fig.[2]3 by explicitly demarcating all regions of phase co¬ 
existence, rather than isolated tie-lines and three-phase 
triangles. Fig. [14] shows all the phase-coexistence re¬ 
gions on a single diagram. We classify the coexistences 
into those with two R phases (Fig. [T3hb two AR phases 
(Fig. [T3bb or R-AR (Fig. [THbb As discussed in [l[, 
each two-phase equilibrium region is associated with a 
metastable promontory, its continuation into a region 
whose equilibrium is three-phase. Hence, the equilibrium 
R-R region of Fig. [T3h is surrounded by a metastable R- 
R region, in which the equilibrium state is R-R-AR. In 
Fig.[T3b. an AR-AR region overlaps two AR-AR-R trian¬ 
gles, which are mutually overlapping (the idea of multiple 
metastable coexistences for a given state point is shown 
in Fig. [5]). All coexistences appearing in Fig. [TSb are 
metastable because AR minima are higher in free energy 
than R. For certain parameter choices, this can be re¬ 
versed [l|. 

R-AR coexistence is available for highly asymmetric 
overall compositions (Fig. [T3b l but has not been the fo- 
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FIG. 13. (Color online) Phase-coexistence regions corresponding to the partial phase diagram shown in Fig. [JJd. Coexistences 
are divided into those containing two R phases, two AR phases, or R-AR coexistence, (a) A narrow two-phase R-R equilibrium 
(dark grey) region abuts two equilibrium three-phase R-R-AR regions (grey). A metastable continuation (red) extends each side 
of the equilibrium R-R region, and example tie-lines are plotted, (b) Mutually overlapping three-phase metastable AR-AR-R 
triangles (reds), and a metastable two-phase AR-AR region (yellow) in which example tie-lines are plotted, (c) Two-phase 
R-AR metastable (blue) and equilibrium (dark blue) regions. The equilibrium region obtains only outside the R-R-AR triangles 
of (a). With the present parameters the equilibrium R-AR regions are close to the edges where numerical evaluation is difficult 
(contrast Ref. 0 ), so they are sketched here by hand (dark blue). 


cus here. Equilibrium R-AR forms a continuation of the 
R-R-AR equilibrium triangles of Fig. [T3k . Due to strong 
segregation for the present parameters, R-R-AR trian¬ 
gles are large; hence, equilibrium R-AR regions exist only 
very close to the edges of the phase diagram. Compare 
[I| in which larger equilibrium R-AR “arms” exist. . 

The finite positive tilt of R-AR tie-lines implies that, 
if each leaflet were imaged separately^ domains of larger- 
than-average appear in both leaflets and be spatially 

colocalised, although the overall leaflet compositions re¬ 
main highly asymmetric [4l|. Thus, for an asymmetric 
bilayer in R-AR coexistence, whether compositional do¬ 
mains in one leaflet “induce” them in the other is related 
to whether the R-AR tie-line is tilted enough for them to 
be detected - a fully horizontal or vertical R-AR tie-line 
would imply compositional domains in only one leaflet. 
This is discussed also in [^ . 

The tilt of R-AR tie-lines is a result of the direct inter¬ 
leaflet coupling B. Without 5, the R and AR free-energy 
minima would be the same height [5T|, and all R-AR 
tie-lines would be vertical or horizontal in (0^, 0^) space. 
Thus, the same direct inter-leaflet coupling that can make 
R phases lower in free energy than AR (rendering R-R 
coexistence equilibrium for symmetrie overall leaflet com¬ 
positions), also gives the ability of one leaflet to “induce” 
domains in the other leaflet where the overall composi¬ 
tions are strongly asymmetric [H, |4l|. Further, R-AR 
tie-line tilt helps explain aspects of the R-R-AR coexis¬ 
tence observed [H . Such observations d, [TgI, |41| therefore 
support a direct inter-leaflet coupling interaction. 

In principle one could prepare a compositionally-pure 
bottom leaflet (for example) and a mixed top leaflet. In 
such a case, the top leaflet can form compositional do- 



FIG. 14. (Color online) Fig. flSk.b and c overlaid to create a 
full phase diagram. 


mains while the bottom leaflet comprises domains of dif¬ 
fering thieknesses (implicitly, tail ordering) despite its 
uniform composition. Via the direct inter-leaflet coupling 
5, the thickness of the bottom leaflet is coupled not only 
to its own composition but to that of the top leaflet. (In¬ 
deed, a one-component leaflet or a one-component hilayer 
can separate into domains of qualitatively differing tail 
order, such as fluid and gel. However, the calculations 
in this paper are explicitly assumed to be far from such 
first-order transitions of the pure lipids - see Footnote 
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[25]). 


The free-energy landscape and available phase equilib¬ 
ria (“phase diagram topology”) can qualitatively change 
for different parameter choices [ol, [l^- For example, 
strong enough B causes AR minima to disappear entirely, 
in which case R-R-AR may still be possible, but AR-AR 
or AR-AR-R coexistence is not. If AR minima exist but 
are of sufficiently different height and curvature to the R 
minima, the AR-AR-R common tangent might disappear 
while AR-AR remains. In Fig. [131 fFe fourfold symmetric 
R-AR arms extend arbitrarily close to the edges of the 
phase diagram, but they could also become truncated 
[l^ such that one pure leaflet (e.g., ^ 0 ) causes the 

bilayer to he outside any coexistence region irrespective 
of the other leaflet’s composition. (In contrast to the 
“induction” of domains described above, this would be a 
case of the pure leaflet’s composition suppressing phase 
separation and thus any domain formation in the bilayer 
as a whole.) The size of the one-phase corners of the 
phase diagram in which no phase separation is possible 
depends quantitatively on the parameters, and they are 
small here due to the strong segregation. 


We believe the phase diagram topology studied here 
(cf. also i) is qualitatively appropriate for experimen¬ 
tal phospholipid systems: R-R-AR coexistence has been 
observed [H , implying that direct coupling should not be 
strong enough to eliminate it; AR-AR coexistence ob¬ 
served in molecular simulation 0 implies that direct 
coupling should not be strong enough to eliminate AR 
minima; we also note that B < J is required to yield pre¬ 
dominant pairwise antiregistration for a laterally homo¬ 
geneous bilayer [l| (as measured in HI), which places an 
upper limit on a realistic value of B that is independent 
of existing estimates of the direct coupling 7 (to which B 
can be mapped Q). Nevertheless, it is true that both ki¬ 
netics and phase diagram topology depend on the values 
of parameters which, at present, can only be crudely esti¬ 
mated in relation to real lipid systems [Il,i0. Refs. 0 
and 0 provide a helpful overview of the broad possibili- 
ties of different phase diagram topologies, albeit without 
considering metastable states and with different (purely 
phenomenological) free energies to that used here. It 
may be useful to use molecular simulation or experiment 
to qualitatively deduce the appropriate phase diagram 
topology by exploring different overall compositions and 
seeing which coexistences are exhibited [ 2 | , perhaps even 
thereby inferring allowed ranges of model parameters. 
This could be complementary to the more bottom-up 
approach of directly estimating the microscopic model 
parameters [l|. 



FIG. 15. (Color online) J = Aa-‘^k^T, B = 0.48 

(0.5, 0.5) overall composition. R and AR structure factors at 

the times indicated (cf. Fig. [6|). 


Appendix B: Registered and antiregistered structure 
factors 


The structure within the top or bottom leaflet alone 
would be measured by 




^rit{h){r,t)exp{iq-r) 


(Bl) 

where angled brackets indicate averaging over indepen¬ 
dent trajectories, and 


Vt(h){r,t) = 


1 if S' at {r,t) in top (bottom) leaflet 

^-1 iiU . 

(B2) 

Averaging over wavevectors of equal magnitude gives 
Such “single-leaflet” structure factors could 
in principle be experimentally obtained by deuterating 
or fluorescently tagging the S lipid of only one leaflet. 
Combining the separate top and bottom structure fac¬ 
tors together yields 

Ssmgie{q,t) = + Sb{q,t)) . (B3) 

5'singie(^: ^) is uot the same as the “overall” structure 
factor of the bilayer S'exp.(^ 5 ^) that would be measured 
experimentally if a given lipid species is tagged the same 
in both leaflets 

2 \ 


^exp. (^5 0 — 


I 


X]^exp.(r,f)exp(fqr • r) 


where 




2 if SS at (r, t) 
1 if SU or US 
0 otherwise . 


(B4) 


(B5) 


We can also define separate registered and antiregis¬ 
tered structure factors Sm{QU) (where m = AR, R) 


Sm{q,t) = 




exp{iq ■ r) 


, (B 6 ) 
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FIG. 16. J = 4a-2teT, B = OASa-^kBT, (0.5, 0.5) overall 
composition. As Fig. [6] but with black lines now showing the 
integrated intensity (zero’th moment) of Sm{qB) instead of 
the characteristic domain lengthscale. 


where 


VR{r,t) = < 


1 if SS at (r, t) 
-1 ifUU 
0 otherwise , 


(B7) 


and 


r]AR{r,t) = < 


1 if SU at (r, t) 
-1 if US 
0 otherwise . 


(B8) 


Writing t/r = |(r^t + ^b) and t^ar = livt - Vh) yields 

Sr + 5'aR = 5'single , (B9) 


i.e., the single leaflet structure factor splits into R and AR 
contributions. Writing ?7exp. = ^r + 1 it can be shown 
that 


5'r(<7, t) = 5'exp.(g, t) . (BIO) 

Hence, the registered and antiregistered structure fac¬ 
tors are in principle accessible via strategic differential 
tagging of the two leaflets. 

From Eqs. 14 and 15, the characteristic lengthscale of 
R and AR domains can be monitored through time. One 
can also view the structure factors directly (Fig. [T5]) to 
compare the position and magnitude of the R and AR 
peaks, or analyse them in other ways. For example, one 
can consider the integrated intensity of Sm{q^ t) (Fig. [16]), 
which is the denominator of Eq. 14. This shows an ini¬ 
tially dominant structural signal from AR domains, fol¬ 
lowed by a downturn which is accompanied by a sudden 
increase in the integrated intensity of SR{qU) as conver¬ 
sion to equilibrium R-R coexistence takes place. 


Appendix C: Linear stability 

As shown in [l|, linear stability analysis of Fq-l 
(E q. 17) yields cj^ax and maximised (over 

wavenumber q) growth rates of the R and AR modes. 



coupling Ja^ / /CbT 


Auj 



FIG. 17. (Golor online) Linear stability analysis as Fig.[4|but 
using V = 0.9 /cbT in the mean-field free energy instead of 
V = 0.6/cbT. 


metastable 



AR-. 


AR-. 


AR-. 







FIG. 18. J = 4a-2/cBT, B = 0.24 (0.3, 0.3) overall 

composition. As Fig. [9l with weaker direct coupling B. The 
bilayer is trapped in AR-AR-R coexistence. 


Fig. m shows Auj = ~ ^mSc- Fig. [17] is as Fig. [4] 

but using V = 0.9 /crT. The stability lines change and 
absolute values of the growth rates are altered, but the 
qualitative landscape of Acj (governing which mode is 
fastest) is not strongly affected. 

The linear stability analysis involves a dimensionless 
parameter the relative mobility for thickness relaxation 
versus lateral diffusion [l|. Figs. 4 and [TTI are calculated 
using ^ = 10 in order to approximate the value Ur = 10 
in simulation. However, in this expected physical regime, 
such that thickness relaxation is much faster than lateral 
diffusion, the results of the stability analysis are only 
weakly dependent on large changes in ^ [l|. 
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metastable 



FIG. 19. J = Aa-‘^k^T, B = 0.24a“^A;BT, (0.5, 0.7) overall 
composition. As Fig. llOi with weaker direct coupling B. The 
bilayer is trapped in AR-AR-R coexistence. 


Appendix D: Trajectories for (0.3, 0.3) and (0.5, 0.7) 
with weaker direct inter-leaflet coupling 

Figs. [18] and [19] show the effect of weaker direct cou¬ 
pling at (0.3, 0.3) and (0.5,0.7) overall compositions re¬ 
spectively. Fig. [18] is as Fig. [9] but with smaller direct 
coupling B. This gives metastably trapped AR-AR-R 
coexistence instead of R-R. Fig. [19] is as Fig. [TO] but 
with smaller 5, and is trapped in metastable AR-AR- 
R instead of R-R-AR. Hence, metastable trapping of AR 
phases can occur even for overall compositions that pro¬ 
hibit perfect (AR-AR) antiregistration. 
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